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The force exerted by a semi-infinite inhomogeneous superconductor with a planar interface to vacuum on a 
magnetic tip is studied theoretically in the absence of external magnetic fields. It is shown that the force has 
a contribution from inhomogeneities due to material defects with unique characteristics. Defects are taken into 
account in the London limit by allowing the mass parameter to vary spatially. The contribution from defects to 
the force is calculated analytically to first order in the deviation of the mass parameter from its constant value for 
the homogeneous superconductor, assuming that the tip is a point dipole perpendicular to the interface, and that 
it does not spontaneously create vortex matter. Random point defects and linear localized defects are considered 
phenomenologically. For each defect type the force dependence on the dipole position coordinates is obtained, 
and the force magnitudes are estimated numerically. The predictions for the dependence of the linear defect force 
on the dipole lateral position are found to agree qualitatively with experiment. 



A magnet and a superconductor placed in the 
vicinity of each other interact through a force 
that results from the changes in the supercon- 
ducting state induced by the magnet. Since in- 
homogeneities in the superconductor due to ma- 
terial defects modify the magnet-superconductor 
interaction, the force must contain contributions 
from defects. This contribution must be present 
in the force on a magnetic tip measured by the 
technique of magnetic force microscopy (MFM) 
PQ applied at zero external field, and can, in prin- 
ciple, be measured if its magnitude is larger than 
the MFM resolution, about lpN. 

For planar superconductor-vacuum interfaces, 
as assumed here, translational symmetry consid- 
erations alone suggest that there is a contribu- 
tion from defects to the force on the tip with 
unique characteristics. For defect-free materials, 
the tip-superconductor interaction energy is in- 
variant under translations parallel to the inter- 
face, and depends only on the tip-superconductor 
distance. This implies that the force on the tip, 
equal to the gradient of the interaction energy, 
is perpendicular to the interface and independent 
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of the tip lateral position coordinates. This is no 
longer true if the superconductor has defects, be- 
cause translational invariance is broken. The tip- 
superconductor interaction energy now depends 
also on the tip lateral position, and so does the 
force, which is no longer perpendicular to the in- 
terface. Thus, the force on the tip has a con- 
tribution from defects with components parallel 
and perpendicular to the interface, both depen- 
dent on the tip lateral position as well as on the 
tip-superconductor distance. This contribution 
can, in principle, be singled out in MFM measure- 
ments since it is responsible for the lateral force 
component and for the dependence of both force 
components on the tip lateral tip position. The 
objective of this paper is to study the properties 
of the contribution from defects to the force and 
to estimate its magnitude using a simple model. 

The general framework for the description of 
inhomogeneous superconductors is the Ginzburg- 
Landau (G-L) model with spatially varying coeffi- 
cients of the G-L free-energy expansion [2] . Here 
the G-L model is used in the London limit, in 
which case only the mass parameter of the G-L 
expansion appears. The effect of the tip mag- 
netic field on the superconductor is to induce a 
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Figure 1. Schematic view of semi- infinite su- 
perconductor and point dipolc M = M z z, at 
r = (Pq,z ). 



screening current and create vortex matter, when 
the tip is close enough to the interface. Here the 
defect contribution to the force is calculated as- 
suming that the tip does not create vortex mat- 
ter. This is the situation of greatest interest be- 
cause, as shown here, the defect contribution re- 
sults solely from to the spatial dependence of the 
mass parameter. 

The model considers a semi-infinite isotropic 
superconductor, separated from the vacuum by 
a planar interface parallel to the x — y plane, 
occupying the region — oo < z < 0. The mag- 
netic tip is approximated by a point dipole per- 
pendicular to the interface, M = M z z, placed at 
r = (x a ,yo, z ) = (p , z > 0) (Fig. [J). The free- 
energy in the London limit, assuming that the 
superconductor is vortex matter free, is given by 
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where r = (p,z), A(r) is the vector potential, 
and A is the penetration depth, related to the G- 
L mass parameter, m, by A 2 = TOc 2 /(47m s e* 2 ), 
where e* = 2e, and n s is the density of Cooper 
pairs. For homogeneous superconductors the 
mass parameter is a constant, mo, and so is 
the penetration depth, Ao- Inhomogeneities are 
described by m dependent on the spatial co- 
ordinates, m(p,z) = to + Sm(p,z). Here 
it is assumed that dm /mo <C 1. The free- 
energy to first order in dm/ mo is written as 
G«(r ) = G (0 \z ) +6G(r ). The first term, 
G^ \ is the homogeneous superconductor free- 



energy, given by Eq. (p} with to = mo, and A 
equal to the equilibrium vector potential in the 
absence of defects, denoted A (r). It depends 
only on zq because of translational invariance par- 
allel to the interface. The second term is the first 
order correction, given by 
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The correction to first order in 5m /mo to the 
free energy, Eq. $2$, is just the change in the ki- 
netic energy of the homogeneous superconductor 
caused by the inhomogeneity, since the screen- 
ing current is given by j sc = — (n s e* 2 /TOoc)Ao. 
The force on the dipole is (r ) = (z Q )z + 
<5F(ro) where the first term is the homogeneous 
superconductor contribution with F^(zq) = 
-dG^{z )/dz and SF(r ) = -V SG{r a ) is the 
contribution from defects. 

The homogeneous superconductor free-energy 
is G(°\zo) = -M z b' z (z )/2, where b' z is the z- 
component of the magnetic field of the dipole 
screening supercurrent at the position of the 
dipole [3]. The calculations carried out here 
use the analytic expressions for b z , and Ao ob- 
tained in Ref.[3]. The force F^°\zo) is repul- 
sive (F^(zo) > 0) and increases as zo decreases, 
typically as shown in Fig. [2l For the dipole at 
Oo = 0,z > 0), A is given by A (p, z; z ) = 
Ag(p, z; z ) 9, and 9 is the unit vector in the di- 
rection of the polar angle. For the dipole located 
at p ^ 0, the vector potential is Ao(p—p , z\ zq), 
because of translational invariance parallel to the 
interface. A typical dependence of Ae on p at 
constant z is shown in Fig. [2l For other values of 
zq and z, Ag decreases smoothly as | z \ increases, 
vanishing for | z |3> Ao, and also decreases with 
increasing zo- 

Next, the contribution from defects to the force 
is calculated for random point defects and linear 
localized defects, assuming in each case that the 
spatial dependence of 5m{p, z) is know. 

Random Point Defects. In this case 5m(r) is 
a random function with zero average satisfying 
< 5m(r)Sm(r') >= 7, n (5(r — r'), where <> de- 
notes average over the randomness [5] ■ According 
to Eq. @ , the force is a random function of the 
dipole lateral position, p 0l with zero average and 
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Figure 2. a) Dependence of (zq) on zq (F^ in 
units of e = (0o/4vrA o ) 2 ) for M Z /4> Q X = 1.0. b) 
Dependence of Ag(p, z; Zq) on p, for po = 0, z — 
— Ao/8, z = 0.3Ao, and M z /</> Ao = 1.0 (Ag in 
units of 4>o/\q). 



root mean square (RMS) fluctuation given by 

<| SF a \ 2 >V 2 = 2^0(^3 ) 1 / 2( ^)2 i?Q(zo ) , (3 ) 
mgA^ O A O 

where a labels the force components (lateral: a = 
p, perpendicular: a — z), eo = (0o/4ttAo) 2 , and 
R a are dimensionless functions, independent of 
M z , given by: 

Rp(zo) = A 2 / d 2 p 1°^ dz | V p al(p, z; z ) | 2 ; 

i? z (zo) = A 2 / d 2 p dz | dal(p, z; z a )/dz | 2 , 
where a , defined by A = (0o/Ao)(M z /</> o Ao)a o , 
is dimensionless. The functions R p (zq) and 
R z (zq) decrease smoothly with increasing z a , 
as shown in Fig. [3^. For a numerical esti- 
mate of <| 5F a | 2 > 1 /2 the value of j m is 
taken from Ref.[2], obtained from collective vor- 
tex pinning theory at low temperatures, namely 
7m/(TOo£ 3 ) ~ 10~ 2 , where £ is the vortex core 
radius. Thus 7m/(m^Ag) - KT 2 (£/A ) 3 . As- 
suming Ao = 150 nm, \q/£, — 10 and zq ~ 
0.3Ao, for which R p , R z ~ 1 (Fig. [3^), gives 
(< (<SF) 2 >)V2 ~ 2 x lO- 2 e o (Af z /0oA o ) 2 
- 2 x 10 -1 (M*/</>oAo) 2 pN, for both the lateral 
and perpendicular force components. Thus, the 
predicted value of (< (SF ) 2 >) 1 / 2 are larger than 
1 pN for M z /4> \ > 2. 

Localized Linear Defect. In this case 8m(r) 
is written as 5m(r) — moXmSoix), where 
Sd{x) — [1 + (a;/a£)) 2 ™] -1 is the defect shape 
function shown in Fig. [3Jd, \m & dimensionless 
constant, ao the defect length scale, and n is a 



positive integer, chosen large in order that the 
boundary between the defect and the homoge- 
neous superconductor is sharp. The assumption 
that 6m/mo 1 is equivalent | Xm 1^ 1- The 
force is given by 

M 

SF(x , zq) = 27re (— -f ) 2 XmL(io, zo) , (4) 

where L = (L X ,L Z ) is a dimensionless vector de- 
fined as: 

L = Ao 2 V / d 2 p j ^ dzal(p, z; z )S D (x+x Q , z). 
Typical results for L x and L z as functions of the 
dipole position coordinates, x$ and zq are shown 
in Fig. |2fc and Fig. [2ji- The highlights of these 
figures are the following. The functions L x and 
L z are non-vanishing only in a range of | xq 
values comparable ajj, and decrease fast with in- 
creasing zq. The former depends on the value 
chosen for ao, and is true in general, except for 
ajj <C Ao. In this case the ranges of L x and L z are 
comparable to Ao, because the range of Ao is Ao, 
but their magnitudes are considerably reduced. 
The function L x oscillates between positive and 
negative values, whereas L z is always negative. 
This behavior results because L x is related to the 
derivative of Sd and L z to Sd itself, and from 
the functional dependence of Sd chosen here. For 
an order of magnitude estimate of the force, the 
same parameters used in the case of random point 
defects are adopted, and Xm ~ 10 _1 is chosen. 
The result is 5F a ~ 27r(Af z /</)oAo) 2 L a pN, where 
a = x, z. The maximum values of the force 
components for zq — 0.3Ao, using the results 
shown in Fig. [3j are 6F X ~ 4O(M z /0 o A o ) 2 piV, 
and SF Z ~ 3 x l0 2 (M z /<j) \ Q ) 2 pN. Thus, for 
zq — 0.3Ao, force values > lpN are predicted 
for modest values of M z /4>q\q, 

The predictions for the force from linear de- 
fects is now compared to MFM force measure- 
ments at zero field reported in Ref.[4]. These 
measurements find that the perpendicular com- 
ponent of the force in the vicinity of a stacking 
fault dislocation in the superconducting material 
depends on the tip lateral position (see Fig. 3 and 
Fig. 4 of Ref.[4]) in a manner similar to that pre- 
dicted here for a linear defect (Fig. (3]i). The in- 
terpretation of these measurements based on the 
results of this paper is that the force behavior is 
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Figure 3. (a) Random defects functions R p (zq) 
and R z {zq) vs. zq- (b) Linear defect shape func- 
tion for old = 0.5Ao, and n = 5 vs. dipole lateral 
position xq- Linear defect functions L x (c), and 
L z (d) vs. xq and zq. 



due to the spatial dependence of mass parame- 
ter caused by the stacking fault dislocation. The 
latter statement is justified by the work reported 
in Ref. [5] showing that defects in the crystal lat- 
tice of the superconducting material cause spa- 
tial dependencies of the G-L parameters, and by 
the fact that the stacking fault dislocation in the 
sample used in Ref. [5] has an approximate linear 
symmetry. The authors of Ref. [4 interpret their 
results differently. They attribute the observed 
behavior to vortices created by unshielded mag- 
netic fields, from the earth and from the exper- 
imental apparatus, pinned by the stacking fault 
dislocation. Both interpretations account qual- 
itatively for the spatial dependence of the force. 
One possible way to distinguish between them ex- 
perimentally would be to revert the tip magnetic 
moment and observe how the force changes. The 
contribution from defects to the force remains un- 
changed, whereas that from vortices changes sign. 

In conclusion then this paper shows that de- 
fects in the superconducting material, as de- 
scribed by spatially dependent G-L mass parame- 



ters, give contributions to the force on the dipole 
with dependencies on the dipole lateral position 
characteristic of the type of defect. The order 
of magnitude estimates carried out here indicate 
that the force from random point defects is prob- 
ably too small to be observed experimentally. On 
the other hand, the force from a linear localized 
defect is estimated to be larger than the resolu- 
tion of the MFM technique, and may already been 
observed. 

The calculations carried out here neglect cre- 
ation of vortex matter by the dipole. This hap- 
pens when the tip is close enough to the inter- 
face, because the dipole magnetic field inside the 
superconductor becomes large enough to create 
vortex half loops [3] . The presence of vortex mat- 
ter changes the defect contribution to the force on 
the dipole. It is expected that, to first order in 
5m /mo, it reduces the magnitude of this contri- 
bution. The reason is that the defect force results 
from the change in the kinetic energy caused by 
the inhomogeneity, and the vortex matter created 
by the dipole generates a current that opposes the 
dipole screening current, reducing the total cur- 
rent and, consequently, the kinetic energy. The 
question is how large this reduction is. To an- 
swer it requires considerations beyond the scope 
of this paper. Work along these lines in under 
way and will be reported elsewhere. 
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